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We present a general approach for evaluating a large variety of three-dimensional Fourier trans- 
forms. The transforms considered include the useful cases of the Coulomb and dipole potentials, 
and include situations where the transforms are singular and involve terms proportional to the Dirac 
delta function (5(r ). Our approach makes use of the Rayleigh expansion of exp(ip- f ) in terms of 
spherical Bessel functions, and we study a number of integrals, including singular integrals, involving 
a power of the independent variable times a spherical Bessel function. We work through several 
fT^ , examples of three-dimensional Fourier transforms using our approach and show how to derive a 

number of identities involving multiple derivatives of 1/r, and Sif). 
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I. INTRODUCTION 



Performing three-dimensional Fourier transforms is a routine task in physics research and in the upper-level physics 
1—1 ■ classroom. A common problem is to transform two-body interaction operators from momentum space to coordinate 
' space and back. Such operators might originate in a momentum-space Feynman diagram analysis but be most 
conveniently evaluated in coordinate space. As useful and representative examples, consider the following transforms 



(~| ' from the text of Berestetskii, Lifshitz, and Pitaevskii 1] in their derivation of the Breit equations Q for electron- 



electron and electron-positron interactions: 



. ■ J _ f <Pp Jp.r Pi /-,, X 

^: ^ J ^-w' 

g : /. = / ^e--^ ^ = Is^.Sif) + ^ (5, - 3x,x,) . (Id) 

The position vector has components f = {xi , 2:2 , X3) = {x,y, z), and Xi = Xi jr. These transforms arise from converting 
the fermion-fermion or fermion-antifermion two-body interaction Hamiltonian from momentum space to coordinate 
space. The first transform shown above is for the usual Coulomb interaction. (Some methods for performing the 
Coulomb transform are shown in Appendix A.) The remaining transforms were evaluated in Ref. |l| by a series of 
^ : tricks. For Eq. (|lbp the method is the simple expedient of writing pi in terms of the derivative of the exponential 
with respect to Xi. 

where di = For Eq. (flcl) they used the exponential derivative to take care of one momentum factor and a 

derivative of with respect to pj to represent the other. The pj derivative was moved to the exponential by an 
integration by parts: 
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The final transform Eq. ()ldp is more subtle. Berestetskii et al. used the derivative trick of Eq. ([2]) twice to write 

/. = -M^. (4) 

They performed the derivatives in Eq. ^ in the usual way for r > to obtain the non-delta function term of Eq. (jldl) . 
They noted that angular averaging amounts to didj — >■ ^6ijd^, which with the point-source Poisson equation 

implies that there must be a term involving the three-dimensional delta function 5{f) on the right hand side of 
Eq. (jld| . They made the implicit assumption that the delta function would only appear in a spherically symmetric 
way, and noted that the angular average of Sij — SxiXj vanishes. With these assumptions, they obtained Eq. (jldp . 
We will have more to say about all of these identities, especially the last, in the remainder of this work. 

We will present a systematic procedure for obtaining transforms of the sort sampled in Eqs. ([1]). We will develop a 
procedure for finding transforms of the form 



(2^ 



(r) - / 7:^e^^">"Pn • • • P., (6) 



when the transform exists. Equivalently we will evaluate the transforms 

InUr) = / ^e^'--p-Yrip). (7) 

The equivalence is due to the fact that any angular function, such as Pi-^^ ■ ■ -Pi^, can be written in terms of spherical 
harmonics Yl^(p). Our method involves expressing the exponential as a Rayleigh expansion ^i5|] 

oo I 

^ip-r ^ ^ ^ (2^ + iy3fXpr)P,{p ■ f), (8) 

e=o m=-t 

where the ii{x) are spherical Bessel functions and the Pt{x) are Legendre polynomials. We use the addition theorem 
for spherical harmonics [6| 

^'^P ■ ^ 27TT ^ Yr{r)Yr{p) (9) 

m— — £ 

to factor the angular dependence of • f) into parts involving the angles of f and p alone. We substitute Eq. ^ 
into Eq. ^ and use the orthonormality of spherical harmonics [7| 

dnpYr{v)Yl!'' {p) = Su,6,nrn' (10) 



to write the transform as 



In,mir) = TT^Yrir) / dpp^+'Mpr). (11) 
The problem has been reduced to the evaluation of the integrals 

/>oo 

i?„,(r) EE / dpp^+'Mpr) (12) 
Jo 

for the relevant values of n and £. 

Transforms of the form shown in Eq. ([G]) can be easily expressed in terms of spatial derivatives, for instance: 

/L-2;.,.....(r) = (-*)X---5,:.y' ^^''''^ = (-')''9,,---d,, (^^y (13) 
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For L > 1 here the derivatives cannot be taken naively near r = and a regulation is required (note the delta function 
in Eq. (|ld|) ). The regulation of derivatives of 1/r has received significant attention recently. |8|-|20| Correspondingly, 
in our direct evaluation of the Fourier transforms the radial integrals Rni{r) don't necessarily converge and must be 
regulated. 

Many of the expressions appearing in this work involve Dirac delta 'functions' or other singular objects. We 
consider these objects to be distributions, or generalized functions, and not functions in the usual sense. [2ll - [25j 
While the values of these quantities might not be well-defined at particular points, their integrals when multiplied by 
an appropriate class of test functions are always well-defined and calculable. A procedure for verifying the correctness 
of formulas involving generalized functions by integration with test functions is described in Appendix B. 

This article is organized as follows. In SectionlUwe describe the evaluation of the integrals i?„^(r). For some values 
of n, t these integrals converge in the strict sense; for some they diverge but useful finite values can be obtained by 
regularization or otherwise; and for others the integrals are proportional to the delta function 5{r). In Section [IIII we 
give results for the Fourier transforms of the forms Inim{r) and In.i-^...i^{r) defined above. In Section |IV] we describe 
several examples and applications, and derive a number of identities involving multiple derivatives of the functions 
1/r, and 5(f). In Appendix VK\ we show how to perform the Fourier transform for the Coulomb potential. Finally, 
in Appendix [B] we confirm a number of our results for generalized functions by integration with test functions. 



II. INTEGRALS INVOLVING SPHERICAL BESSEL FUNCTIONS 



We need values for the integrals Rni{r) defined in Eq. p^ . For convergent integrals we can change the variable of 
integration to obtain Rniir) = xm/r"'^^ where 



Xnt= / dxX-^+'jtix). (14) 
JO 

In order to know what values of n and i lead to convergent integrals, and in order to get values, we will need to look 
carefully at the spherical Bessel functions, [l^ll^] The basic definition of ji{x) is in terms of the usual Bessel function 
(of the first kind): 

= J-^Ji+i/2{x), (15) 
V 2x 



which is actually defined for complex values of £. A power series expansion for ji{x) is 
We will need only integral values £ = 0, 1, • • • , in which case 



j,{x) ^i-xY [-—] . (17) 

- ax J X 



The first few spherical Bessel functions are: 

jo{x) = , (18a) 

X 

, , sin a; cos a; „, , 

3r(x) = — , (18b) 

x'^ X 

. . / 3 1 \ cos X , , 
^^^""^ " xj '''''' ^^""^ 

We see that they are relatively simple functions with oscillatory behavior that fall off for large x. The small x behavior 
is less easy to see from the explicit forms, but from the power series expression in Eq. (IT6| it follows that 

Mx)-j^ + Oix^-^) (19) 



for small x. The large x asymptotic behavior is given by 



, , sin(x — n£/2) 
je{x) « ^ L^. (20) 
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FIG. 1: Values of n and £ for which Rni{r) = dpp"^^ji{pr) is defined. In the region on the left, with —{i + 3) < n < — 1 
(indicated by small dots), R„i{r) and Xn« = r^^'^ Rniif) = dx x"'^^je{x) are actually convergent. The solid dots show 
integral values of n, I in this region. In the region on the right, with —2 <n < i (indicated by dashed lines), Rni{r) and Xnt 
can be given values either by extending the domain of Eq. (|2ip . by reciprocity Eq. (|24|l . or equivalently, by introduction of a 
cut-off as in Eq. (|27|l . Circles show integral values of n, £ in this region. On the line £ = n, Rne(r) is proportional to the delta 
function 5(r) as in Eq. (|32|) . Integral values on this line are represented by squares. 



It follows that in order for the integral Xnt to exist, it must be that —1 < n + 2-\- ^ lest there be a divergence at small 
values of x, and also n < — 1 to avoid a divergence at large values of x. We note that only the first and third among 
the transforms of Eqs. ([T]) lead to convergent radial integrals. The tables give the value [28[ 

X„. = 2"+iV^^L^ (21) 

for the integral ([T4|) as long as n is in the range of convergence: —(£ + 3) < n < —1 (for the case with n and £ real). 

Transforms (jlbl) and (jldl) have n ~ —1 and n = respectively, which are outside the range of convergence for 
our radial integral, so we will have to extend that range. This extension can be done both formally and physically- 
fortunately, the various approaches are consistent with one another. As a formal method of extension we can simply 
allow Eq. (I^T|) to define the value of Xni for values of n and £ for which it would otherwise not be defined. This is 
the same approach that is used to define values for the divergent integrals of dimensional regularization. (29l - [33[ This 
procedure allows us to extend the range oi n to —{£ + 3) < n < i; when n = £ the value given by the formula Eq. ((2T|) 
vanishes and hides the important fact that Ru{r) contains a delta function singularity at r = 0. The values of n and 
£ for which integral Rne{r) can be defined in this way are represented in Fig. [T] 

Another formal approach to the extension of the allowed range of n makes use of the completeness relations for 
spherical Bessel functions [33 - [36{ 

2r^ f°° 

(5(r-a) = / dpp'^je{pr)ji{pa). (22) 



TT 







An interesting derivation of Eq. (j22p based on the Rayleigh expansion has been given by Ugincius. [s?] Starting from 
Eq. (p2)) we divide by a"+^ and integrate over a: 

1 da . , . 2r2 p da 1 , , , 

Changing variables now to z = pr, x — pa we obtain the reciprocity relation 



dzz^+^Mz) / dxx-("+i)j,(a;) (24) 



TT 

2 ' 

Relation ([M)) is of a formal nature as the range of convergence — (£ + 3) < n < — 1 of the z integral only overlaps with 
the range —2 < n < of the x integral in the region —2 < n < —1. Nevertheless, if we use Eq. (PT|) to perform the x 
integral in its allowed range, we find for the z integral the result 

dzz-+'Mz) = ^ 2-+'V^-^-jl-L = xn,, (25) 

2x-(„+3),£ r (^) 
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which agrees with Eq. (|2ip . This allows us to extend the range of n for which Xni is defined to the whole region 
-(^ + 3) < n < £. 

As a physical regularization we make use of the expectation that the infinite extent of the range of integration is 
an idealization and will, in one way or another, be cut off for large values. We model this cut-off by a mathematically 
expedient exponential decrease, and generalize Rn(.{r) to [s^ 

dpe ^Pp^'+^Mpr) = ^,+1 T{£+l) ^ ^^^^ 

for positive A. In the limit A the integral (1^51) becomes 



A^O+./n r"+3 r(^) -^"+3 ^ ^ 



consistent with Eq. ([?!]) but convergent for all — (i' + 3) < n. 

We will also need useful values for the integrals Rni{r) of Eq. (|12p with n — £. As a simple first example we work 
out the case with £ ^ Q. Starting from Eq. ([2^ in the limit a — > 0, wc find 

27-2 r°° 

Kr) = / dpp^joipr). (28) 



7^ JO 

This result has been noted in Ref. [39l It is a consequence of the usual Fourier representation for the delta function 

= / jSf^-"' = 3^ i" /' = i /°° ""^^ = i ""'"fe-' <^'" 

(where u = cosOp as in Appendix A) after use of the delta function identity S{f) — S{r)/ {Anr'^) (see Appendix B for 
a proof). In order to obtain the needed generalization, we rewrite Eq. (f22|) remembering that the integral vanishes 
unless r ~ a: 

(5(r-a) = ^ dpp^jt{pr)jt{pa) (30) 

and take the a -> limit using Eq. (HH): [i^l 

or 

o/^ r, i+2-> . ^(2^+1)!!., ^ 27r2(2£ + l)!!^_^ 

Ru{r)^J^ dpp'+^jpXpr) ^ \^t+2 '^(^) ^ (32) 

We now bolster this formal derivation with a physical one based on the inclusion of a cut-off. The cut-off version of 
Ru{r) is the integral 

(We have included a pre-factor for convenience.) Then Ri{X; r) is a representation of the delta function. It is properly 
normalized: 

poo 

drRi{\;r) = l (34) 

for all positive values of A. jHI The functions i?^(A; r) have the general shape shown in Fig. O [i^l For each value of 
t, Ri{X; r) has a single maximum of width 0(A), height 0(1/A), and position a distance of 0(A) to the right of r = 0. 
The functions Ri{X;r) satisfy the sifting property 



lim / dr/(r)i?,(A;r) = /(0) (35) 



for all continuous functions /(r) that vanish sufficiently rapidly at large r, and so form representations of S{r) for 
each value of £ (including non-integral values of £ as long as —3/2 < £). We notice that Xne = when n ^ £ from 
Eq. (mi) is consistent with xu = i?«(l) oc 6(1) = from Eq. (15^ . 
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RfiX,r) 




FIG. 2: The functions Re{X,r) that form representations of the delta function S{r), shown for £ — 0, £ = 3, and £ = 10 with 
A = 0.04. These functions all have unit area and a single maximum whose height is 0(1/A), width is 0(A), and position is a 
distance of 0(A) to the right of the origin. 



III. RESULTS FOR THE FOURIER TRANSFORMS 



We can now present results for the three-dimensional Fourier transforms considered in this paper. The transforms 
of the form of Eq. (O that we can do are 

InUr) = J ^e^'^-'p-YriP) = ^^rrir), (36) 

when — (£ + 3) < n < i, and 

- J ^''''p' - = /^^^<5(f )rr(f), (37) 

which is the form appropriate when n — £. Results for these transforms have a simple form because each contains 
only a single value of angular momentum £, and so involves only a single spherical Bessel function. 

A more general result for transforms like that of Eq. p6p can also be found that generalizes the number of spatial 
dimensions to N and allows for arbitrary complex values of n. [i^ The general result has a form much like that of 
Eq. (|36p (but generalized to N dimensions) except for the exceptional lines £ = n — 2k (where k = 0, 1,2, . . .) on 
which the transform contains the delta function or its derivatives, and £ — —n ^ N — 2k (where fc = 0, 1, 2, . . . ) on 
which the transform has a somewhat complicated form containing In r and powers of r but no delta functions. 

The transforms © arc more complicated than those of (O because they involve multiple values of £. We reduce 
Eq. ([6]) to a sum of terms like Eq. ([7]) by expanding each angular factor Pi^ ■ ■ ■ Pi^ in terms of spherical harmonics as 

1 or f 

pn---p^. = J2 J2 ct-rjr{p). (38) 

t=L m=-t 

where £ in the sum takes the values L, L — 2, . . . , 1 or (depending on whether L is odd or even). This expansion can 
be understood by remembering that each factor of has £ = 1, as follows from the fact that each component of p can 
be expressed linearly in terms of (p) , and the addition law of angular momenta implies that the maximum value 
of £ obtained from L factors of unit angular momentum is L. Parity tells us that only odd or only even values of £ 
contribute-as p^^ • • ■ pij^ has parity (—1)^ and Y"''{p) has parity (—1)^, and all terms of Eq. must have the same 
parity. Explicit values of the Cf™.j^ coefficients can be found using the orthogonality of the spherical harmonics: 

= J Yr ip) Pn--- P^. • (39) 
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We turn now to a study of the angular momentum decomposition of Pi^ ■ ■ ■ Pi^ ■ From our previous discussion it is 
ear that Pi^ ■ ■ ■ Pi^ contains contributions of angular mo] 
to be that part oi Pi^ ■ ■ -pi^ of angular momentum i.e. 



clear that Pi,^ ■ ■ ■ Pi^ contains contributions of angular momenta L, L — 2, etc., down to 1 or 0. We define (pi-^ ■ ■ ■ pij^ )g 



{p,,---p,,)^^ ^ cf™.,,rr(p). (40) 

m— — £ 

An explicit expression for {pi^ ■ ■ -Pi^)^ can be obtained by using Eq. ([5^ in Eq. (PU)) and the addition theorem for 
spherical harmonics: 

{Pn■■■P^.)i^ I dn,,p'^^---p[^ Yr{p')Yr{p)^{2i+l) I ^p',^---pi^Pe{p'-p). (41) 

m——£ 

The method is to write out Pe{p' ■ p) as a polynomial of order t and perform the angular integral using [4^ 



An 



(jV +7)!! ■ ■ ■ '^'"-1'" + P™)(JV-1)!! terms ' (^2) 



(The integral vanishes for odd iV, and the subscript indicates that there are a total of {N — 1)!! terms in the sum 
over permutations.) Using Eqs. (HT|) and it is not difficult to obtain (p^^ • • -pi^)^ for £ = and £ — I because 
the corresponding Legendre polynomials are so simple. (The first two Legendre polynomials are Po{p' ■ p) ~ 1 and 
Pi{p' ' P) = p' ' P = P'jPj with an understood sum over j from 1 to 3.) We find that 

(P'H ■ ■ ■ P^Jo = J^Pn-- -pI = (^^+7)1! ('^'1^2 ■ ■ ■ K-i^L + perms) ^^^^^^ , (43a) 



[Pii ■■■PiL)i = 3 / -j^ p,^ ■ ■ ■ Pi^pj pj = 3 (p,, ■ ■ ■ p,^Pj)q Pj 

3(5L,odd 



(L + 2)!! (^''^ ('^^^^^^ ••■'^^--1^- + P™)(L-2)!!terms+Pe™s)^t„,^, ■ (^Sb) 



This direct approach becomes increasingly cumbersome for higher values of i. However, the results obtained in 
Eqs. (|43p can serve as initial values for an inductive proof of a general expression for (p^^ • • •Pij^)i . [45j 

The most convenient way to obtain {pi^ ' ' ' PiL)e^ j least for small L, is by straightforward construction. We 
explain the method here and give several examples. We note first that the quantity {pi-^ ■ ■ -pi^)^ is symmetric in all 
indices. Furthermore, the part with maximal angular momentum, {pi^ ■ ■ -pij^)^, is traceless: 

iP^l■■■P^L)L^^L-l,^L={P^l■■■P^L-2)l ^=0, (44) 

because an object of angular momentum L cannot be constructed from a combination of only L — 2 objects each 
having angular momentum one. The condition of tracelessness allows us to construct {pi-^ ■ ■ ■ Pi^)^ relatively easily 
for each particular value of L. The L ~ 2 result is immediate: 



2 1 

{PtiPt2)2 = PiiP^2 - 3^n»2- (45) 



The remainder is the £ = contribution: 



so that 



2 1 

(p»iPi2)o = 3^n»2> (46) 



PilPi^ = {P^lP^2)2 + [PilMl ■ (47) 

For L = 3 we start with Pi^Pi2Pi3 and add a symmetric term with two fewer momentum unit vectors but the same 
three indices with an initially undetermined coefficient: 

(^11^2^43)3 ^PtlPl2Pt3 + A{p^^5^2^^ +Pl2Sl3ti +P43^n«2) ■ (48) 



Contraction over 12 and yields 



13/3 "1223 

SO that A = —1/5. So for i = 3 we find 



= {PhPi2Pt3)3 ^^2i3 =Pii+ A (5p,J , (49) 



3 1 

{PilPi2Pi3)3 = PilPi2Pi3 - g iPilSi2i3 +Pi2^i3il +^^3^*1^2) (50) 



as the maximal angular momentum term, and the remainder is the £ = 1 contribution: 



3 1 

{PilPi2Pi3)l = g {Pil^l2i3 +Pi2^l3il +Pi3^lll2) ■ (51) 



By construction, we have the angular momentum decomposition: 



PllPt2Pl3 = {PilPl2Pl3)l + {PtlPl2Pl3)\ ■ (52) 



The term with four momentum unit vectors can be decomposed in an analogous way: 



{PtiPi2Pt3Pu)t = PtiPi2Pt3Pt4 - y {PtiPi2^t3i4 + perms)g t^j.^^ + {S^,^2S^^u + perms)3 ^^^.^^ . (53) 



The £ ^ 2 term of 

PiiPi2Pi3Pii comes from subtracting the appropriate traces from the parts of PiiPi2Pi^Pii 
{PiiPi2Pi3Pii)\ that are quadratic in momentum unit vectors: 



The £ = Q part is what remains: 



{PtlPt2Pt3Pti)t = y ( iP^lP^2)l ^^3^4 + Pe™s)g j^^^^. (54) 



(p»iP*2P*3P»4)o = ^ {^^l^2Si3u + perms)3 j^^, . (55) 



As a final example, we give the results for L = 5 

g {P^^P^2P^3K^, + PCrms)io terms + ^3 



{Pi4Pi2Pi3PiAPi^)l = P2iPz2Pz3PuPi5 " „ (^^1^2 '^*4*5 + pcrms)^^ ^^^^^ + — {Pz, Si^^J,,i, + pcrms)^^ j^^j^^ ,(56a) 



{PhPi2P^3PuP^, )l = li iP^lPi2P^3)l '^u»5 + perms) (56b) 

{P^^P^2P^3P^4P^r^\ = ^ (Pn ^»2»3'5»4»5 + Pems) (56c) 

We note that angular decomposition works for vectors with non-unit length just as well as for unit vectors. One 
simply has that 



(Pn • ■ ■ )£ = P^ iPn ■ ■ ■ PiL )e ■ (57) 



So, for instance. 



where 



and 



PilPi2Pl3 = {PllPi2Pi3)l + {PllPt2Pl3)l (58) 

3 

{PtlPl2Pl3)3 ^PllPl2Pt3 - Y (Pn^l2t3 +P^2^^3^1 +P23<5n*2) (59) 

3 P 

{PllPi2Pl3)l = Y (Pil'^»2»3 +Pt2Sl3ll +Pl3^ilt2) ■ (60) 

We are now ready to obtain results for the transforms Inyi-^^...i^{f ) of Eq. ([6]). First, for the component {pi-^ ■ ■ -pi^)^ 
of angular momentum £, one has 

re^ P (P^l■■■P^L)i ^ 7r-^-;n^ ix^l■■■x^L)e ^ (61) 
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TABLE I: Summary of the three-dimensional Fourier transforms evaluated in this work. The functions <l?(p) and 'I'(r ) are a 
transform pair as defined in Eqs. (|65p . In each of these transform pairs one can make the replacement Y[^{p) — >■ (p^^ ■ ■ ■ Pij^)f 
and Yi!^{r) (x^j ■ ■ ■Xi^)f to obtain a different representation of the angular part of that transform. The value of n is 
restricted to the range —{£ + 3) < n < £. The value of Xn.e is given in Eq. (|21l) . 

<I>(p) ^(f) 

/vrip) /m^5{r)Yr{f) 

^_iY(2i±m^2nfS{p)Yr{p) r'Yrir) 



valid when —{£ + 3) < n < £. The result (pT|) comes from Eq. ([55]) . the expression (PH)) giving {pi-^ ■ ■ -pi^)^ in terms 
of Yf^{p), and an analogous expresssion for {xi^ ■ ■ ■ Xi^)^ in terms of Y^™(f): 

(x,,-..£,J,^^ ^ Cl^..,Yr{f). (62) 

When n = £ the transform contains a delta function: 



(27r)3 

by use of Eq. ([57)) . Transforms In-^i^-.-i^ir) of the function 

-p"5](Ph---^J,^ (64) 
I 

are done through use of superposition. 

The Fourier transforms that we have done are summarized in Table HI which shows the transform pairs $(p) and 
^(r ). The functions $(p) and ^'(r ) are related by 



fp 
(2^ 



4>(p) = J d^re-'P-^^{f). (65b) 

(The radial part of this transform can be expressed as a Hankel transform, (46l . |47| so many additional transform pairs 
are also known.) 

IV. APPLICATIONS AND CONSEQUENCES 

As first examples of the use of our Fourier transform formulas we return to the transforms of Eqs. ([T]). The first 
three are immediate applications of Eq. (ICTI) . One has for la the result 



/.= //^e-l = /,^ = -^, (66) 
J (27r)3 p2 2tt^ r inr 

as X-2.0 = 7''/2 and the angular decomposition of L momentum vectors when L = is ()[] = 1. Terms 7^ and Ic can 
be done just as easily. For Id we find 

= J {{PrPjfo + {PrP3)l] = 5{^){XrXj)l + ^^ix,Xj)l = ^'5y'5(f ) - (^X^Xj - ^6^j^ (67) 
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because xq.2 = 37r/2. According to Eq. Q we can write the transform Id as a double spatial derivative of 1/r, which 
leads to 



^^^A = - + ^ - l^ij ] ■ (68) 



We note that Eq. is the sum of two terms: one that comes from the usual derivative formulas and that holds 
when r ^ 0, with the other contributing only when r — 0. This identity was emphasized by Frahm, Q and has 
many applications. Its immediate applications are for the calculation of the electric and magnetic fields for electric 
and magnetic dipoles. [22, 48, 49] The electric and magnetic dipole potentials are given in terms of the electric and 
magnetic dipole moments p and rfi as Ve = p ■ rjr^ and Am = fh x r/r^ . The corresponding fields are 

„ / „ 1\ 47r , , 3(p ■ f)f — r^v 
E = -VVe = ipjd,-j - -yp(5(r ) + ^ (69) 

and 

B = V X A„ - V X TO X V- = -wSI- + V TO • V - = — TO (5(f) + ^ '- , (70) 

\ r ) r \ / r 3 r° 

where V is the gradient vector (Vi = di) and is the unit vector along the i*^ coordinate axis. (We have used 
the vector identity Ax (^B x = B (^A- — C (^A- B^ .) The electric and magnetic dipole fields have identical 

structure outside the source but are quite different at the position of the dipole. The electric dipole is produced by 
separated charges of opposite sign, with the dipole moment p pointing from the negative charge to the positive one. 
The field between the charges points from the positive charge to the negative one as represented by the negative 
sign on the pS{r) term in Eq. (j69|) . The magnetic dipole is produced by an infinitesimal current loop with all of the 
flux of the dipole flowing up through the loop in the direction of the dipole moment to, corresponding to a positive 
sign on the mSif) term in Eq. (|70l) . The delta function term in Eq. ([70]) is the origin of the hyperfine splitting for 
atomic S states. [10, [HI In particular, in the ground state of atomic hydrogen it leads to the 21cto line (53 - l54| that 
is important in radio astronomy. 

We now turn to a number of additional Fourier transforms that give interesting results. First we consider the 
transform of ^ (pj^ ■ ■ ■ Pj,-)'l- On the one hand, this transform can be written in terms of derivatives: 

(03-"-'^ (Pn ■ --P^.t = i-^f (% • • • dj.t -tr^ (71) 
while on the other hand, we can evaluate the transform using Eq. (j6ip : 



(27r)3^ p2 ■■■Pjk)k- 27^2 '^^■J'l '"^ik)k^ y'^) 

where from Eq. pip we obtain Xk-2.k ~ ^(2fc — 1)!!- We find then for the derivatives of l/r the formula 

(% • ■ ■ 5.J^ \ = ^''^'^^''r'^" i-n ■ ■ ■ ■ (73) 

We now write this out explicitly for a few low values of k. For k = the identity is trivial (l/r = l/r) where we recall 
that Oq = 1 and (—1)!! = 1. For fc = 1 we obtain the expected result: 

a,i^-^. (74) 
r 

For fc = 2 wc obtain 

which, with use of the Poisson equation ([5|) , leads again to Eq. ([68| . The fc = 3 result can be expressed as 

didjdk^ = {xiXjXk)l - ^{di6{r)5jk + djS{r)6ki + dk6{r)Sij), (76) 
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an identity also obtained by Frahm. Q We point out that although the general k*'^ derivative of 1/r contains delta 
functions, the particular combination {dj-^ ■ ■ ■ dj^)i^ ^ with maximal angular momentum has none. 

The second new transform we consider is that of ^ {pj-^ ■ ■ ■ Pjk)^- This transform can be represented in terms of 
derivatives of l/(27r^r^) (the transform of 1/p), or evaluated explicitly using Eq. (IFTI) : 



Making use of Xk-i,k — S'^fc!, we find that 

1 (-l)''2*=A:l k 
i^n ■ ■ ■ )k^= (%i ■ ■ • )fc • (78) 

Again, the combination of derivatives of 1 /r^ having maximal angular momentum gives a result without delta func- 
tions. Relation ([75]) is trivial for k = 0. For fc = 1 it gives 

d^--'-^^ (79) 

and for fc = 2 

- ls,.,d'^ ^ = ^ (id, - i<5.,) . (80) 

These formulas are easy to verify for r > 0; the new result is that no extra delta functions appear. 

The final new transform we consider is that of {pj-^ ■ ■ -pjj.)^. As before, this transform can be represented either in 
terms of derivatives (this time of S{f)), or explicitly by use of Eq. We find that 



fp_ .p.w .,.(2fc + l)!! 

(2^)3 



(-0'(% •••%): -5(0= / ^e^P-(p,,•••p,J:=^^•^-^<5(r)(%...^,J^ (81) 



The new delta function identities are 



(% • • • d,jl m = ^ (%. ■ • ■ J' m- (82) 



and 



For fc = 1 and fc = 2 these are 

d,S{r) = -—5{r) (83) 
r 

- \5,,d'^ S{r) = ^ (^x,x, - S{r). (84) 

These can be checked using the method of Appendix B. 

Appendix A: Fourier transform of the Coulomb interaction 

We first consider the Fourier transform of the Coulomb potential from coordinate space to momentum space. Direct 
evaluation in spherical coordinates doing first the (trivial) (j) integral and next the 6 integral leads to 



/ 



d^re-iP-^ = -/ drr d6l sin 6*6- 'P'™"*' = - / drrl ■ =-/ drrsm{pr). (Al) 



Anr 2 7o 7o Jo \ ipr J P Jo 



(It is convenient to use the substitution u = cos9 with limits — 1 < u < 1 for the 9 integral.) This integral has a 
problem at large distances: it doesn't converge. We repair this 'infrared' divergence by the insertion of a convergence 
factor in the original integrand: 1/r — )• e^^^ /r with A a small positive number that will be taken to zero at the end 
of the calculation. With this alteration, the transform (IA1[) becomes 

d3^g-.pr ^^_L Tdr (e-(^-'f - e-(-^+'^P>) - = ^1—. (A2) 



Airr 2ip Jq ^ I 2ip \X — ip X + ip J p^ + 
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On taking the A — ^ limit this shows that the transform of l/(47rr) is The inverse transform dial) then takes 

l/p^ to l/(47rr). 

The original momentum space to coordinate space transform of Eq. ([Taj) can be performed as a contour integral. 
We retain the momentum space version of the infrared convergence factor and evaluate the (j)p and Op integrals: 

"^'^ re^P-^^^^ ^-^ r dp( + ie^^- - e'^^^) . (A3) 



(27r)3 p2_^A2 Sn^r J„ \p - iX p + iX 

We note that the integrand is even in p so we can extend the integration range to — oo < p < (X) if we divide the result 
by 2. We use the residue theorem to evaluation the integral, closing the e*'"' term in the upper half plane and the 
g-ipr ^Qj-i^ |;jjg lower half plane. The two contributions are identical: the final result is e~^^ /(Anr), as expected. 
The convergence factor was convenient to move the poles off of the real p axis. Without the convergence factor the 
integral (|A3p would be 



(27r)^ p^ 27r^ Jq pr Airr 

a familiar (and tabulated) integral obtained by the one-dimensional Fourier transform of the rectangular pulse function 



Appendix B: Verification of identities involving generalized functions 

In this appendix we describe our means of confirming identities involving generalized functions and illustrate with 
proofs of several equalities from the text. Suppose A and B are two generalized functions and consider the assertion 
A = B. We adopt the 'physicists' proof approach of Frahm [8]] and say that A ^ B holds if 

(frF{r)A = j (frF{r)B (Bl) 

for all test functions F{r) that are smooth (infinitely differentiable) and that decrease sufficiently rapidly at large 
distances. The generalized functions considered in this work have singularities at only a single point, which (by 
translational invariance) we have taken to be the origin. So an identity like A = B can be tested at points different 
from the origin by traditional means. We let the integration region implicit in Eq. (IBip include the origin, and break 
it up into two parts: a small sphere (of radius R) surrounding the origin, and an outer part. As long as A — B 
away from the origin the integrals over the 'outer part' will be automatically identical and we needn't consider them 
further. So the identity A = B holds when A = B away from the origin and 

/ dPrF{f)A^ / SrF{r)B (B2) 

J B(R) Jb{R) 

in the limit i? — )■ where the integration region is the ball B{R) of radius R surrounding the origin. The integral over 
B{R) can be written as 

(fir = [ drr^ [ dQ, (B3) 

B{R) Jo J 

where dfl is the element of solid angle, and by convention we always perform the angular integrals first. For some 
of the identities to be considered, one of the generalized functions is a gradient. Suppose then that A = Q^A. We 
evaluate the integral over B{R) according to 

d^rF{f)dkA ^ [ d^rdk{F{f)A)~ [ d^r{dkF{f))A 

B{R) JB{R) J B{R) 



(f dVLR^XkF{f)K- [ d^r{dkF{f))A, (B4) 

Js(R) Jb(R.) 



IS{R.) JB(R.) 

where the integral in the first term is over the sphere S{R) of radius R that is the surface of B{R). The surface 
element of S{R) is dS — dSf = dilR^r. The integral identity used in transforming the first term in Eq. (jB4p is a 
generalization of Gauss's Theorem. |55j] The test functions F{f) are infinitely differentiable and so can be expanded 
in Taylor series around the point r = 0: 

F{r) =Fo+ XaFa + ^XaXbFab + ' ' ' (B5) 
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where Fq = F(0), Fa = daF{0), Fab = dadbF{0), etc. It follows that dkF{r) = Fk + XaFak + \xaXbFabk + ■■■■ 
In Sec. [IT] it was asserted that the three-dimensional delta function can be written as 

The proof of this statement is now straightforward. The left hand side (LHS) of Eq. (IB1|) is simply F(0) by use of 
the sifting property of the three-dimensional delta function. The right hand side (RHS) of Eq. (jBl|) is 

RHS = — / dr5{r) [ dQ F(r,e.(f>) = F(0), (B7) 
47r 7o J 

since the one-dimensional delta function has unit area in the region < ?' < i? for any positive R and F{0, 9, (f) — F(0) 
for any 9, (p. The equality (jB6[) is thus confirmed. 

As a second example, we consider the proof of Eq. ([80)1 . In this case the left hand side of Eq. ([80t is a gradient 
9fcAfe where 

Afe = (s.kdj - ^Sijdk) \ = — % (SikXj - ^SijXk] ■ (B8) 



3 

The first term on the LHS of Eq. dHTJ (the first term of Eq. (IB4|) is 



LHSi ^ (j^ dn R^XkF{r) ^ j (^S.^Xj - ^(5yXfe j ^ ^ j> dVl {f^ + XaFa + ]^XaXbFab + • • • ^ {x^Xj)l . (B9) 

Now the {xiXj)^2 term has angular momentum 1 — 2 and will vanish when multiplied by _Fo (which has ^ = 0) or 
Xa (which has £ = 1) and averaged over angles. The XaXbFab term is the first in the Taylor expansion to survive 
the angular average, but it contributes at order R and so disappears in the i? — )■ limit. Higher order terms in 
the expansion of F{r ) are even higher order in R and so do not contribute either. The second term on the LHS of 
Eq. ((BT|) (the second term of Eq. ((B4)) ') is 



LHSs = - drr^ J i^k + rXaFau + • • ■ ) {-^^ {5^kXj - ^5ijXt^ 



(BIO) 



The angular average eliminates the product of Fi (which has £ = 0) with Xj or x^ (which both have I = \). The next 
term in the series for F{r) contributes at order R and so vanishes as do all higher terms in the expansion. The RHS 
of Eq. dm]) is 

RHS = J^ drr^ J dnF{r )^{x^Xj)l. (Bll) 

The first term in the expansion for F{f) that survives the angular averaging is •^XaXbFab, which leads to a negligible 
contribution of order R. Identity (IMl) holds for r > by explicit calculation, and it holds near r = because 
LHS = RHS (both being 0), so Eq. (iSOl) is an identity among generalized functions. 

As a final example, we consider the delta function identity given by Eq. (I84p . The Aj, quantity for this case. 



Afc = yS,kdj - -Sijdkj 5{f) = l^SikXj - -SijXkj S{f), (B12) 

itself contains a delta function and so vanishes unless r = 0. It follows that the surface term, LHSi, is zero. The 
second LHS term is 

LHS2 = - J d^'r {dkF{f)) = - j d\ {Fk + XaFak + • • • ) (^5^kXJ - '^S.^xt^ 5{r) = Fj - ^^^Fkh- (B13) 

A similar calculation shows that the RHS is also Fij — ^SijFkk, completing the proof. 
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